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ABSTRACT 

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with 
different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. 
Each grid covers masses ranging from 5 to 60 Mq and initial rotation rates between and about 600km/s. 

To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective 
overshooting by using the observed drop in rotation rates for stars with surface gravities log^ < 3.2 to determine the width of the 
main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in 
the Large Magellanic cloud. 

We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the 
occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the 
most massive fast rotators. 

We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and 
the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that 
the consideration of surface abundances can resolve this degeneracy. 

Key words. Stars: abundances - Stars: evolution - Stars: early-type - Stars: rotation - Stars: massive 
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I 1 . Introduction tion jMitalas & Falkll 19841: [Gonzalez et al.ll2005h . However, two 

major issues have plagued evolutio nary models of massive stars 
O Massive stars can be considered as cosmic engines. With their ^^^^j j^^^y. ^-^-^^ ^^^^ 1^^^ (IChiosi & Maedeii^986l) . We 
— ' high luminosities, strong stellar winds and violent deaths they concentrate here on the role of mixing in massive stars as, on 
drive the evolution of galaxies throughout the history of the ^^-^ sequence, the elfects of mass loss remain Umited in the 
universe. Even before galaxies formed, massive stars are be- considered mass and metalhcity range, 
lieved to have played an ess ential role in re-ionizing the uni- 
verse (iHaiman & Loebl[T997l) . with important consequences for Concerning thermally driven mixing processes, the occur- 
its subsequent evolution. Massiv e stars are visible out to large r^nce of semiconvection and the related question of the appro- 
distances in the nearby universe (iKudritzki et al.ll2008h . and as P^^te choice of the criterion for convection do not alter the main 
ensembles even in star-forming galaxies at high redshift (e.g., sequence evolution of massive stars significantly (|Langer et al.| 
lDouglasetal.ll2(i9h . Their explosions as supernovae or gamma- \12M)- However, the efficiency of convective overshooting in 
ray bursts shine through a major fraction of the universe, probing massive main-sequence stars is still not well known. In partic- 
intervening structure s as well as m a ssive star evolution at the "lar, the cool edge of the main-sequence band is not well deter- 
lowest metallicities (ISavagliol 12005 lOhkubo et al.l iSl) . It is mined, as stars in observational samples are found in the gap pre- 
therefore of paramount importance for many areas within astro- dieted between main-sequence stars and blue core-heliurn burn- 
physics to obtain accurate models of massive star evolution. In ers (e.g. [Fitzpatrick & Garmanv||1990t [Evans et al.||2006|). Here 
this paper we present an extensive grid of evolutionary models we follow a promising approach to determine the cool edge of 
of rotating massive stars focusing on the main-sequence stage, *e main sequence band, and thus the overshooting efficiency 
which is at the root of understanding their subsequent evolution, from the rotational properties of B-type stars (see [ Hunter et al. | 

Basic concepts about massive main-sequence stars are | 2008b| ; |Vink et alj | 2010|) . 

long established and provide a good guide to their ex- Rotationally-induced mixing processes have been invoked 

pected properties, most importantly the mass-luminosity rela- in the past decade to explain the surface enrichment of 
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some massive main-sequence stars with the products from 
hydrogen burning, in particular nitrogen (iHeger et al.l l2000t 
[Meynet & Maeder 2000; Maeder 2000). Rotationally-induced 
mixing was further used to explain the ratio between O-stars 
and vari ous types of Wolf-R ayet stars at different metalHci- 
" I2005h . the variety of core collapse su- 



ties jMey net & Maeder 
pemovae (IGeorgv et al. 200 9) and the evolution towards long 
gamma-ray bursts (lYoon &'Lan ger 2005). However, in view of 
the fact that all the mentioned observed phenomena have alter- 
native explanations, a direct test of rotational mixing in massive 
stars appears of paramount importance. 

In view of the recent r esults from the FLAMES Survey of 
Massive Stars dEvans et al.ll2005i) . which provided spectroscopic 
data for large samples of massive main sequence stars, we pur- 
sued the strategy of computing dense model grids of rotating 
massive main sequence stars, and predict in detail and compre- 
hensively their observable properties. This evolutionary model 
data 0, which was calibrated using results from the FLAMES 
Survey, is tailored to provide the first direct and quantitative test 
of rotational mixing in rapidly rotating stars dHunter et al.l2008aL 
[2009; Brott et al. 2010, from here on Paper II). 

Our new model grids deliver predictions for many observ- 
ables, with a wide range of potential applications. Besides he- 
Uum and the CNO elements, we show how the light elements 
lithium, beryllium and boron are affected by rotational mixing, 
as well as the elements fluorine and sodium. We also provide 
isochrones which predict complex structures near the turn-off 
point of young star clusters. In Sec. |2] we describe the set-up 
and calibration of our massive star models. Sec. |3]andHprovide 
the evolutionary tracks in the HR diagram and the correspond- 
ing isochrones, and the surface abundances of all considered el- 
ements. We end with a brief summary in Sec. |5] 



2. Stellar evolution code 

We use a one-dimensional hydrodynamic stellar evolution code 
that takes into account the physics of rotation, magnetic fields 
and mass-lo s s. Th is code has been described extensively by 
Heger et all (200dl). Re cent improvement s are presented in 
Petrovic et al.l ( l2005h and lYoon et al](l2006l) . 

To allow for deviations from spherical symmetry due to rota- 
tion in this one-dimensional code, we consider the stellar prop- 
erties and structure equatio ns on mass s hells that cotTespond 
to isobars. According to e.g. IZahnI (1 19921) turbulence efficiently 
erases gr adients along isobaric su rfaces and enforces shellular 
rotation (lMevnet&Maedej|1997l) allowing us to use the one- 
dimensional approximation. The effect of the centrifugal accel- 
er ation on the stellar structure equ ation s is considered according 
to lKippenhahn & ThomasI (Il970h and lEndal & Sofial (Il976h as 
described in lHeger et al 1(120001) . 



2. 1 . Transport of chemicals and angular momentum 

All mixing processes are treated as diffusive processes. 
Convection is modeled using the Ledoux criterion, adopting a 
mixing-length parameter of aMLT = 1.5 (lBohm-Vitenselll958t 
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1 (iLanget 



lLangedll99ll) . Semi-convection is treated as in 
(Il983h adopting an efficiency parameter ffsEM 
Il99lh . We take into account convective core-overshooting us- 
ing an overshooting parameter of 0.335 pressure scale heights 
i.e. the radius of the convective core is equal to the radius given 



' Our model data is made available through VizieR at 
|http://cdsarc. u-strasbg.fr/cgi-bin/. . .i 



by the Ledoux criterium plus an extension equal to 0.335 Hp, 
where Hp is the pressure scale height evaluated at the formal 
boundary of the convective core (see Sec. l2.4l for the calibration). 
Furthermore, we consider various instabilities induced by rota- 
tion that result in mixing: Eddington-Sweet circulation, dynam- 
ical and secular s hear instability, a nd the Goldreich-Schubert- 
Fricke instability (Heger et al."2000(). 

Transport of angul ar momentum is a lso tr eated as a diffusive 

Jrocess following lEndal& Sofial (119781) an d iPinsonneault et aP 
1989h as described in iHeger et al.l (l2000l) . The turbulent vis- 
cosity is determined as the sum of the diffusion coefficients for 
convection, semi-convection, and those resulting from rotation- 
ally induced instabilities. In addition, we take into account the 
transport of angular m omentum by magnetic fields due to the 
S pruit-Tayler dynamo (ISpruitll2002l) . implemented as described 
in iPetrovic et al.l ( l2005l) . We do not consider possible transport 
of chemical elements as a result of the Spruit-Tayler dynamo 
because its validity is still controversial (iSpruit2006i) . We note 
that, while the dynamo process itself was c onfirmed through nu- 
merical calculations bv lBraithwaitd (l2006l) . it is still unclear un- 
der which co nditions it can op erate and whether it is active in 
massive stars (IZahn et al .1120071) . Heuristically, angular momen- 
tum transport through magnetic fields produced by the Spruit- 
Tayler dynamo, as the so far only mechanism, appears to repro- 
duce the rotation rates of white dwarfs and neutron stars quite 

well (Heger et al. 2005; Suiis et al. 20081). 

In agreement with IMaeder & Mevne3 (l2005l) we find that 
magnetic fields keep the star near rigid rotation throughout its 
main sequence evolution, suppressing the mixing induced by 
shear between neighboring layers. The dominant rotationally in- 
duced mixing process in our models is the Eddington-Sweet cir- 
culation, a large scale meridional current that is caused by ther- 
mal imbalanc e in rotating stars between pole and equator (e.g. 
lTassoull[T978l) . 

Some of the diffusion coefficients describing rotational mix- 
ing are based on order of magnitude esti mates of the relevant 
time and length scales (iHeger et al.ll2000l ). To consider the un- 
certainties efficiency factors are introduced, which need to be 
calibrated against observational data. The contribution of the ro- 
tationally induced instabilities to the total diffusion coefficient is 
reduced by a factor fc - 0.0228 (see Sec. l2.5l for the calibration) 
while their f ull value enters in the expression for the turbulent 
viscosity (see lHeger et al.l2000i for details). The inhibiting effect 
of chemical gradients on the efficiency of ro tational mixing pro - 
cesses is regulated by th e parameter /„, see iHeger et al.l JToOOl) . 
We adopt = 0.1 after lYoon et al.l (l2006l) who cahbrated this 
parameter to match observed surface helium abundances. 



2.2. Mass loss 

We updated the treatment of m ass loss by stellar winds by im- 
plementing the pr escription of Vink et al] (l2000l 1200 1 1) , based 
on the method by Ide Koter et al. I (119971) . for winds from early 
O- and B-type stars. This mass loss recipe predicts a fast in- 
crease of the mass-loss rate as one moves to lower tempera- 
tures near 22 000 K. This increase is related to the recombina- 
tion of Fe IV to Fe m at the sonic point and is commonly re- 
ferred to as the bi-stability jump. These rates are derived for 
12.5 kK ^ Teff $ 50 kK. We note that the formulae on 
the hot and cool side of the jump have not been derived for 
the intermediate temperature range between 22.5 and 27.5 kK. 
In this range, we thus perform a linear interpolation. In or- 
der to accommodate for a strong mass-loss increase when ap- 
proaching the HD limit, we switch to the empirical mass loss 
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rate of lNieuwenhuiizen &de.Taged(ll990h. when the Vink et al 
rate be comes smaller than that from lNieuwenhuiizen & de JageJ 
(II99OI) at any temperature lower than the critical temperature for 
the bi-stability jump. This ensures a smooth transition between 
the two mass loss prescriptions. This strategy also naturally ac- 
counts for the increased mass loss at the second bi-stability jump 
at -12.5 KK. 

To account for the effects of s urface enric hment on the stellar 
winds, we f ollow the approach of Yoon et al. (2006). The mass- 



loss rate of IVink et al.l (l2000l l200lir is employed for stars with 
a surface helium mass fractio n, Fs. below 0.4. We in terpolate 
between the mass-loss rate of IVink et al.' ('2000', 2W^ and the 
Wolf-Rayet m ass-loss rate of | Hamann et al. (X995} reduced by 
a factor of 10 dYoon et al.ll2006h for 0.4 <~r, < 0.7. When the 
helium surface abundance exceeds Ys > 0.7, we adopt the Wolf- 
Rayet mass-loss rate. 

Rotational mass lo ss enhancement is included according to 
lYoon&Lan^ (120051) . 

2.3. The initial chemical composition 

We adopt three different initial compositions that are suitable for 
comparison with OB stars in the Small and Large Magellanic 
Cloud (SMC, LMC) and a mixture that is tailored to the Galactic 
sample of the FLAMES survey of massive stars (lEvans et al.l 
I2005h . to which we will refer in this paper as the Galactic (GAL) 
mixture for brevity. In contrast with several previous studies we 
do not simply adopt Solar-scaled abundance ratios. Even though 
such mixtures may be sufficient to study the overall effects of 
metallicity, they are not accurate enough for dir ect comparison 
with o bserved surface abundances. For example. iKurt & Dufouj 
(Il998h find that the carbon to nitrogen ratio in HIT regions in 
the LMC and SMC is considerably larger than the solar ratio. 
In the stellar interior, carbon is converted into nitro gen, which 
can rea ch the surface due to mixing and/or mass loss. lBrott et al.l 
(l2008h found that the surface nitrogen abundance can be en- 
hanced by a factor of 11 with a realistic initial SMC mixture. A 
Solar-scaled mixture with the same iron abundance would lead 
to enhancement of only a factor of 4, for the same initial param- 
eters. This demonstrates the need for tailored initial chemical 
mixtures. 



iHunter et al.ll2007l l2009h . Also for the SMC these authors find 
no evidence for any systematic differences. The Fe abundance 
for the SMC is taken from a study of A-supergiants by Mem] 
(119991). For all the ot her elements we adopt solar abundances by 
lAsplund et"al1 (|2005|) scaled down by 0.4 dex for the LMC mix- 
ture and 0.7 dex for the SMC mixture. The initial abundances 
are summarized in Tab. [1] These mixtures have metallicities of 
Z = 0.0047 for the LMC and Z = 0.0021 for the SMC, where 
Z is defined as the mass fraction of all elements heavier than 
helium (see Table|2|. 

Choosing a chemical composition suitable for comparison 
with stars in the Milky Way is not straight forward due to the 
chemical gradients within the disk of our galaxy. We adopt 
an initial Mg and Si abundan ce based on t he ana lysis of the 
Galactic cluster NGC6611 by iHunter et al.1 (l2007l). For C we 
adopted the NLTE corrected value from Hunter et al.l (l2007l) . N 
and O abund ances are based on the average of Hll-regions, see 
IHunter et al.l (r2008b„ .2009i) and references therein. The Fe abun- 
dance is taken from measurements of A-supergiants by IVennI 
([T995I) . F or all other element s we have adopted solar abun- 
dances by lAsplund et alj ( l2005h . The resulting metallicity of our 
Galactic mixture is Z = 0.0088, w hich is lower than the so- 
lar metallicity of Z = 0.012 found by lAsplund et al.l (l2005h and 
the solar neighborho od metallicity of Z = 0.014 measured by 
iPrzvbilla et al.l (l2008l) . We point out we assumed the stellar wind 
mass loss rates as well as the opacities in our models to depend 
on the iron abundance, for which the spread in measurements is 
much smaller than for the total metallicity. Our wind mass loss 
rate scales with (Fesmf /Fep)"'^^, where we use the iron abun- 
dance from lGrevesse et all (1 19961) as solar reference, for consis- 
tency reasons with older models. 

For helium we assume that the mass fraction scales lin- 
early with Z between the primordial helium mass fraction of 
Y = 0.2477 (Peimbert et al. 2007b at Z = and Y = 0.28 at 
the solar value (Grevesse et al. 19961). We adopt the OPAL opac- 
ity tables dlglesias & Rogers,, 1996i) using (Fesurf /Fcq) x Zq to 
interpolate between tables of differe nt metallicities. The solar 
reference values are again taken from lGrevesse et aT l (IT996l . 



Table 2. Resulting hydrogen (X), helium (Y) and metal (Z) mass 
fractions for the chemical mixtures used in our models. 



Table 1. Initial abundances for C, N, O, Mg, Si, Fe adopted in 
our chemical compositions for the Magellanic Clouds and the 
Galaxy (see Sec. 12. 3b . 





C 


N 





Mg 


Si 


Fe 


LMC 


7.75 


6.90 


8.35 


7.05 


7.20 


7.05 


SMC 


7.37 


6.50 


7.98 


6.72 


6.80 


6.78 


GAL 


8.13 


7.64 


8.55 


7.32 


7.41 


7.40 



Notes. All other elements are solar abundances dAsplund et al.ll2005l) 
scaled down by 0.4 dex for the LMC and by 0.7 dex for the SMC. 



For the initial abundances of C, N and O we use data from 
HII regions by Kurt & Dufour ( 1998). Surface abundances for 
these elements in B-type stars may no longer reflect the initial 
composition due to nuclear processing and mixing. For Mg and 
Si we use measurements of B-type stars as these stars are not 
expected to significantly alter the abundances of these elements 
during their main-sequence evolution. The B-type stars in dif- 
ferent fields and clusters in the LMC sho w only small varia- 
tions in Mg and Si of less than 0.02 dex dTrundle et al.l 120071: 





X 


Y 


Z 


LMC 


0.7391 


0.2562 


0.0047 


SMC 


0.7464 


0.2515 


0.0021 


GAL 


0.7274 


0.2638 


0.0088 



2.4. Calibration of overshooting 

Mixing beyond the convective stellar core can occur when con- 
vective cells penetrate into the radiative region, when they "over- 
shoot" the boundary due to their non-zero velocity. This effect is 
accounted for through a parameter a which measures the exten- 
sion of the affected region in units of the local pressure scale 
height. In general, overshooting leads to larger stellar cores re- 
sulting in higher luminosities and larger stellar radii, especially 
towards the end of the main sequence evolution. 

Various attempts have been undertaken to con strain the value 
of this parameter Using eclipsing binary stars, ISchroder et al.l 
(119971) derived values between 0.25 and 0.32 for stars in the 
mass range of 2.5 to 7 Mq. iRibas et alj ( |2000|) and lClaretl (|2007|) 
found a in the range of 0.1 to 0.6, with a systematic increase of 
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Fig. 1. Projected rotational velocity versus surface gravity for 
stars in the FLAMES survey near LMC clusters NGC 2004 
(triangles) and Nil (circles). Color coding indicates the stel- 
lar mass. We use the sudden transition at logg - 3.2 (vertical 
line) to calibrate the amount of overshooting of our 16 model 
(black line), see Sec. 12.41 For comparison we plot several other 
models with different initial masses and rotation rates (dotted 
lines). 



the am ount of overshooting with the stellar mass. iBriquet et al.l 
(l2007h used asteroseismological measurements and obtained 
a - 0.44 + 0.07 for the /?-Cephei d Ophiuchi. Another way 
to constrain the overshooting parameter is to fit the width of 
the main-sequence band f or clusters with tum-ofi" mas ses be- 
low approximately ISMq dMermilliod & Maedei^'1986'). Using 
non -rotating evolutiona r y mod els, Mermilliod & Maeder (198a) 
and iMaeder & Mevnetl (Il987l) found an overshooting param- 
et er between 0.25 and .3. This is supported by findings 
of iNapiwotzki et aTl (Il991h . In general, however, this method 
does not work well for masses above 5...IOM0, since the red 
edge of the main sequence band remains mostly unidentifiable 
(IVink et al.ll2()T0h . 

In this work we re-calibrate the overshooting parameter for 
our models using data of the FLAMES survey of massive stars. 
We use the fact that the properties of stars near the end of their 
main-sequence evolution depend on their core size. For larger 
overshooting parameters, the end of the main sequence shifts to 
lower effective temperatures and lower surface gravities. After 
the end of their main sequence evolution, the rapid stellar ex- 
pansion results in a strong spin down of the envelope. This phase 
of evolution occurs on the Kelvin-Helmholtz timescale, which is 
very short compared to the nuclear timescale. Observing stars 
during this phase is unlikely. 

In Fig. [U we plot the projected rotational velocity against the 
surface gravity for sta rs in the LMC samp le of the FLAMES 
survey (adapted from iHunter et al.ll2008"bl) . Stars with surface 
gravities of logg > 3.2 show a wide range of projected ro- 
tational velocities, usini, in clear contrast with the stars with 
lower surface gravities which have projected velocities of about 
50km/s or lower. We interpret this transition as the division be- 
tween main-sequence stars that are bom with a wide range of 
rotational velocities and a diff'erent, more evolved, population 
that has experienced significant spin down. The surface gravity 
at which this transition occurs seems to be independent of mass, 
at least within the range available in this sample. To show this 
we color-coded the masses of the observed stars in Fig. \T] 
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Fig. 2. Similar to Fig.[T] but now for the SMC (top) and Galactic 
(bottom) samples of the FLAMES survey. The typical error in 
logg is about 0.15 dex. The drop in rotation rate, which may 
be interpreted as the end of the main sequence, occurs around 
logg = 3.2 (vertical line) independently of metallicity. 



In an alternate scenario discussed b\ jVink et al.l (|20I0|) . the 
population of slow rotators might be main-sequence stars, that 
have spun down due to increased mass loss (bi-stability break- 
ing). However, this would require a huge overshooting, in order 
to extend the lifetime of main sequence stars at temperatures 
below the bi-stability temperature of about 22 kK. With the cur- 
rent overshooting calibration, bi-stability breaking plays a role 
only well above 30 M©. The nature of the evolved slow rota- 
tors in Fig.[T]remains puzzling. The fact that all of them have a 
strong nitrogen surface enhancement might propose that they are 
post-red supergiants which acquired the nitr ogen enhancement 
through the first dredge-up (IVink et al.ll20I0l) . but they might as 
well be products of binary interaction or of so far unidentified 
physical processes. 

For the overshooting calibration we use a stellar model with 
a mass and rotation rate that are representative for the entire 
sample. Th e typical mass is about I6M0 (based on evolutionary 
masses; see lPaper H) and the mean projected velocity of the stars 
in the sample with logg > 3.2dex is (ysin/) - llOkm/s. This 
corresponds to an average rotational velocity of (v) = 142 km/s 
assuming a random orientation of the spin axes, i.e. (sin /) = tt/A. 
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Throughout the main-sequence evolution the equatorial velocity 
is expected to be nearly constant due to the reduced effect of 
stellar winds in the low metallicity environment of the LMC and 
due to coupling of the expanding envelope to the contracting 
core. Therefore, we use the typical velocity as an estimate of the 
average initial rotational velocity. 

The full black line in Fig. [1] shows a model with the typi- 
cal mass and rotation rate. The overshooting parameter has been 
adjusted to a = 0.335, such that the end of the main sequence 
coincides with the drop in rotational velocities at \ogg - 3.2. 
Fig. [U demonstrates that the adopted amount of overshooting, 
derived for a 16 Mq star, appears to be valid for the mass range 
10 to 20 Mq. However, it is unconstrained outside this range, 
as for example the 30 Mq tracks in Fig. [T] have no observed 
counterparts near the TAMS. We therefore use an overshooting 
parameter of a = 0.335 for our entire grid. Given the uncertainty 
in the surface gravities (~0.15 dex), we estimate the uncertainty 
of the overshooting parameter in the mass range 10 to 20 M© to 
be 0.1. 

A priori we do not expect the overshooting parameter to 
depend on metallicity. To check this we show the projected 
rotational velocities against surface gravity for the SMC and 
Galactic sample of the FLAMES survey in Fig. [2] The drop in ro- 
tation rate occurs at the same surface gravity within the measure- 
ment errors. The Galactic sample shows a second feature around 
log? = 3.5 caused mainly by stars below 5Mo. Duftonetaj] 
(120061) already described this stars as over-luminous for their 
age and a general disappointing agreement of the clusters HR- 
diagram with isochrones. When we vary the initial composition 
from SMC to the Galactic mixture, we find that the end of the 
main sequence for our 16 M© model shifts from log^ ~ 3.3 
to 3.1, adopting the same value for the overshooting param- 
eter. This shift is also within the measurement uncertainties. 
Therefore, adopting one fixed value for the overshooting param- 
eter in our entire model grid is a reasonable assumption. 
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Fig. 3. Nitrogen surface abundance versus projected rotational 
velocity for stars in the LMC sample of the FLAMES survey 
dHunter et al.l2009l) in the Nil (blue) and the NGC 2004 (green) 
fields. Stars with evidence for a binary companion from radial 
velocity variations are plotted as triangles. Black solid lines indi- 
cate re gions that are helpful as a reference in later discussion(see 
also Paper Ili) . Evolutionary tracks of 13 M© models are plotted 
for various rotational velocities, shifted by a factor ;r/4 to ac- 
count for random inclination angles. The red line marks the ni- 
trogen abundance reached at the end of the main sequence for 
13 Mq models adopting a mixing efliciency of = 0.0228 (our 
calibration). Thin dotted lines show the location of the line, if 
different values for fc are adopted, e.g. fc - 0.015 (lower line) 
and fc = 0.03 (upper line). 



2.5. Calibration of tine rotational mixing efficiency 

As mentioned in the beginning of Sec.|2]we treat all mixing pro- 
cesses as diffusive processes, where the diffusion coefficient is 
taken as the sum of the individual diffusion coefficients for the 
various mixing processes. The contributions of all rotationally 
induced mixing processes are reduced by a factor, , before they 
are a dded to the total diffusion c oefficient for transport of chem- 
icals. iPinsonneault et al] (Il989h was the first to introduce this 
parameter In order to explain the solar lithium abundance they 
needed to reduce the efficiency of rotationally induced mixing 
process es by a factor fc - 0.046 . Based on theoretical consid- 
erations IChabover & ZahnI (1 19921) proposed fc ^ 0.033, which 
was adopted in the models of Heger et all (l2000h . We note that 
we cannot compare our efficiency parameter directly with these 
values due to differences in the implementation and the consid- 
ered specific set of mixing and angular momentum transport pro- 
cesses. 

We calibrate the mixing efficiency for our models using the 
B-stars in the LMC sample of the FLAMES survey. We aim to 
reproduce the main trend of the observed nitrogen surface abun- 
dances with the projected rotational velocity for the nitrogen en- 
riched fast rotators. In Fig. |3]we present the data, indicating dif- 
ferent regions or box es which hav e been numbered for reference 
(iHunter et al.ll2008al and lPaper iJ) . The majority of the sample, 
about 53% of the stars, is not significantly enriched in nitrogen, 
given the typical error of 0.2 dex for the nitrogen abundance 
measurements. This group is indicated as Box 5 in Fig. |3] Fast 



rotators with significantly enhanced surface abundances (> 7.2 
dex; see Box 3) constitute 18% of the sample. A detailed discus- 
sion of the stars in Box 1 and 2, which deviate from the main 
trend, and of observational biases an d other fa cto rs that might 
influenc e this numbers are provided in lPaper lH and lHunter et al.] 
(I2008ah . 

For the calibration we use stellar models with a mass that 
is typical for the stars for which we have nitrogen surface abun- 
dance measurements, i.e. 13 M©. We adjust the mixing efficiency 
such that the nitrogen surface abundances reached at the end 
of the main-sequence evolution follow the trend of the stars in 
Box 3 while we make sure that the slow rotators do not en- 
rich more than 0.2 dex to match the stars in Box 5. We find 
that this is obtained for a mixing efficiency of /c=0.0228. In 
Fig.|3]we plot evolutionary tracks for several 13 M© models with 
different initial rotation rates. Adopting this mixing efficiency, 
a model with the typical mass and average rotational velocity 
(fzAMS = 142 km/s) reaches a surface nitrogen enhancement of 
7.2 dex at the end of its main sequence. To illustrate the sensi- 
tivity of the nitrogen surface abundances to the adopted mixing 
efficiency, we depict the surface abundances reached at the end 
of the main sequence for fc = 0.015 and fc = 0.03 (dotted lines 
and grey band). As we have no reason to believe that the mixing 
efficiency fc depends on metallicity, we adopt the same value for 
the SMC and the Galactic grid. 
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Fig. 4. Surface abundance of nitrogen (bottom row) and helium (top row) reached at the end of the main sequence (colors and 
contours) in our models with different initial masses, rotational velocities and compositions. Each evolutionary sequence of our 
grids is represented by one dot in this diagram. See Sec. l2.6l for a description of the initial parameter space of our grid and Sec.|4] 
for a discussion of the abundances. 



2.6. Grid of evolutionary tracks and isocinrones: initial 
parameters and description of provided data 

Using the calibration and the three initial compositional mix- 
tures described above we computed a grid of evolutionary mod- 
els. We adopt initial masses between 5 and 60 and initial 
equatorial rotational velocities (wzams) ranging from to about 
600km/s. For each initial composition we computed about 200 
evolutionary tracks. The spacing between the initial parameters 
was chosen to provide a proper coverage of the parameter space 
and to provide extra resolution for massive fast rotating stars 
whose evolution depends sensitively on the initial rotation rate 
(see Fig.HJi. 

The results of our calculations are made available onUne 
through VizieR in table format. For each evolutionary track we 
provide the basic stellar parameters at different ages such as 
the mass, effective temperature, luminosity, radius and mass-loss 
rate as well as the surface abundances of various elements. We 
also list the central and surface mass fractions of the most im- 
portant isotopes of these elements. In Table[3]we summarized all 
the provided data. The time steps are chosen to resolve the evo- 
lutionary changes of the main stellar parameters. Tracks with 
higher time resolution are available from the authors upon re- 
quest. 

Based on these evolutionary tracks we generated sets of 
isochrone s, using our population synthesis code STARMAKER 
(iPaper iJ) . The isochrones cover the evolution until central hy- 
drogen exhaustion for stellar masses between 5 and 60 M©. We 
provide them for ages between and 35Myr at intervals of 
0.2 Myr, for initial rotational velocities between and 540 km/ s 
with steps of 10 km/ s and for three different initial compositions: 
SMC, LMC and a Galactic mixture. We provide the basic stellar 
parameters and surface abundances as listed in Tab. |4] 



Table 3. Description of the online tables containing the results 
of our evolutionary calculations as a function of time, for given 
initial composition, initial mass and initial rotation rate. 



stellar 


surf. 


surface & core 


parameters 


abun 


isotope mass fractions 


Age (t [yrs]) 


H 


'H 






Mass (m [Mq ]) 


He 


3 He 


■*He 




Eff. temp. (Teff[K]) 


Li 


'Li 






Luminosity (log(L/Lo)) 


Be 


'Be 






Radius (R [Rq]) 


B 


tog 


"B 




Mass loss rate (M [ MJyv]) 


C 








Surface gravity (log g) 


N 


14n 






Surface velocity (usmf [km/s ]) 





16q 


"O 


18q 


Rotation period (P [days]) 


F 


19p 






Critical velocity (Ucrit [km/s ]) 


Ne 


20Ne 


-'Ne 


22Ne 


Eddington factor (Fe) 


Na 


^^Na 








Mg 


^^Mg 


25Mg 


26Mg 




Al 


26 Al 


2'Al 






Si 


28Si 


'-"Si 


»Si 




Fe 


56pg 







Notes. The critical velocity is computed as f^^t" GM/R(l - Fe), where 
Fe denotes the Eddington factor due to electron scattering. 



3. Evolutionary tracks and isochrones 

3.1. Evolutionary tracks 

Rotation affe cts evolutionary tracks in the Hertzsprung-Russell 
diagram (e.g. llVIaeder & ]VIevnelll206oi and references therein). 
The centrifugal acceleration reduces the effective gravity result- 
ing in cooler and slightly less luminous stars. However, rotation 
also induces internal mixing processes which can have the op- 
posite effect, leading to more luminous and hotter stars. Which 
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Table 4. Description of the online tables containing stellar pa- 
rameters along isochrones computed for different initial rota- 
tional ve locities as provided by our population synthesis code 
dPaper iJ) . See also Tab. [3] 



stellar 


surf. 


parameters 


abun 


Age (t [yrs]) 


H 


Mass (m [Mo ]) 


He 


Effective temperature (Teff [K]) 


B 


Luminosity (log(L/Lo)) 


C 


Radius (R [Rq]) 


N 


Mass loss rate (M [MJyi]) 


O 


Surface gravity (log g) 


Ne 


Surface velocity (Vg„yf [km/s ]) 


Na 


Rotation period (P [days]) 


Mg 


Critical velocity (ucrii [km/s ]) 


Si 


Eddington factor (T^) 


' H (mass fraction) 


initial mass and velocity 


''He (mass fraction) 



of these effects dominates depends on the initial mass, rotation 
rate, metallicity and the evolutionary stage. In this section we 
describe the effects of rotation on our models. 

In rotating stars, the radiative energ y flux depends on the lo- 
cal effective gravity (Ivon Zeipelll 19241) and temperature and lu- 
minosity become latitude dependent. The resulting thermal im- 
balance drives large scale meridional currents. In most stars, 
rotating at moderate rotation rates, the expected latitude de- 
pende nce of temperature and luminosity is week (iHunter et al.l 
l2009l) . In models close to critical rotation, luminosity and tem- 
perature differences are expected to be more significant and may 
even give rise to polar winds a nd equatorial outflows caused by 
critical rotation (iMaederll 19991) . Temperatures, luminosities and 
gravities given for the models presented in this paper are surface 
averaged values. 

Mass dependence of rotational mixing: For intermediate mass 
main-sequence stars the main effect of rotation is a reduction 
of the effective gravity by centrifugal acceleration. Fig. |5]shows 
that, during the main-sequence evolution, the tracks of the fast 
rotating 5 Mq stars are about 2000 K cooler compared to the non- 
rotating counterparts. Rotational mixing does affect the surface 
abundances of those elements that are so fragile that the rela- 
tively low temperatures in the stellar envelope are sufficient to 
induce nuclear reactions on them. However, the mixing is not 
efficient enough to significantly alter the structure of the star. 
As a consequence, the effect of rotation on the corresponding 
evolutionary tracks in the Hertzsprung-Russell diagram remains 
limited. 

In rapid rotators more massive than about 15 M©, the effect 
of rotational mixing becomes more important than that of the 
reduced effective gravity: towards the cool edge of the main se- 
quence band, the rotating stars become more luminous than the 
nonrotating ones. The main effect of rotation is an effective in- 
crease of the size of the stellar core, similar to the effect of over- 
shooting (Sec. 12.4b . This can be seen in Fig. |5] when one com- 
pares the tracks of the 15 and 16 Mq models near the end of their 
main-sequence evolution. The larger core mass in the fast rotat- 
ing models results in a higher luminosity. 

In the most massive stars at low metallicity (> 16 for 
the SMC, M > 19 Mq for the LMC), and for the extreme rota- 
tors shown in Fig. |5] mixing induced by rotation is so efficient 



SMC 

6.5 I ; ; , 




50 40 30 20 10 

Teff [kK] 



Fig. 5. Evolutionary tracks for various masses and metallicities. 
Red tracks are the non-rotating models and in blue we show 
models rotating initially at about 550 km/s. For claiity, we alter- 
nate different line styles and scaled the temperature ranges such 
that the maximum of each plot is visible. 



that almost all the helium produced in the center is transported 
throughout the entire envelope of the star. During their main- 
sequence evolution they become brighter and hotter, evolving 
up- and bluewards in the Hertzsprung-Russell diagram. After all 
hydrogen in the center has been converted to helium, the star 
contracts, reaching effective temperatures of up to 100,000 K. 
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6.2 
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Teff [kK] 

Fig. 6. Evolutionary tracks for 15, 25 and 50 Mq for SMC, LMC 
and Galactic composition (from top to bottom), for a range of 
initial rotational velocities (see legend). 



This type of evolution i s referred to a s "quasi-chemically ho- 
mogeneous evolution' ' (iMaeded fT987t lYoon & Langeil l2005t 
IWooslev & Hegeill2006h . The blueward and redward evolution 
after the exhaustion of hydrogen in the center is analogous to the 
"hook" at the end of the main sequence seen in tracks of non- 
rotating stars. Note that all models have been computed beyond 
central hydrogen exhaustion, in most cases up to the ignition of 
helium. This allows us to later compare the main sequence sur- 



face abundances with those resulting from the first dredge-up in 
the red supergiant regime. 

Metal I icily dependence of rotational mixing: Fig. |6] depicts 
evolutionary tracks of models with various initial rotation rates. 
The most striking feature in these diagrams is the bifurcation 
of the evolutionary tracks occurring at high masses and low 
metallicity, most clearly visible at SMC metallicity. Stars that 
rotate faster than a certain threshold are so efficiently mixed 
that they evolve almost chemically homogeneously. Stars that 
rotate slower than this threshold build a chemical gradient at the 
boundary between the convective core and the radiative enve- 
lope. This gradient itself has an inhibiting effect on the mix- 
ing processes, strongly reducing the transport of material from 
the core to the envelope. The minimum rotation rate required 
for chemi cally hom ogeneous evolution decreases with increas- 
ing mass (lYoon et al. 2006). 

At high metallicity, rotational mixing is less efficient. In ad- 
dition, mass and angular momentum loss due to stellar winds 
becomes important, slowing down the rotation rate and there- 
fore the efficiency of the mixing processes. The fastest rotating 
stars at high metallicity initially evolve blue- and upward in the 
Hertzsprung-Russell diagram, see in the lower panel of Fig. |6] 
However, the combined effects of spin down by a stellar wind 
and the build up of an internal chemical gradient reduces the 
efficiency of internal mixing processes. The star switches onto 
a redward evolutionary track, similar to that of a non rotating 
star. However, its larger core mass results in a higher luminosity 
compared to the slower rotating counterparts. 

3.2. Isochrones 

While for a given metallicity, a classical isochrone can be rep- 
resented by a single line in the Hertzsprung-Russell diagram, 
the isochrones of rotating stars span an area for a given age and 
initial composition. This is shown in Fig.|7] The isochrone con- 
structed from nonrotating evolutionary models is plotted in red. 
The isochrone corresponding to the fastest initial rotational ve- 
locity for which the models do not yet follow a blueward evo- 
lution in the HRD is plotted in black. The isochrone with the 
slowest initial velocity that shows a clear bluewards evolution 
is shown in green. In blue and yellow we have selected two 
isochrones from the transition region that may help the reader 
to asses the sudden transition from classical to chemically ho- 
mogeneous evolution. 

For rotational velocities above 350km/s, when the most 
massive stars undergo chemically homogeneous evolution, the 
isochrones deviate strongly from the nonrotating case. The max- 
imum spread in effective temperature occurs around 4 Myr (sec- 
ond row in Fig. |7]i. At lower metallicity the isochrones split into 
more clearly separated branches. At SMC metallicity (e.g. the 
third panel in the first column in Fig. IT) the isochrone based on 
models rotating initially at 500km/s moves straight to the blue. 
In contrast, the comparable LMC isochrone returns to the red for 
stars above ~5OM0. This behavior is directly related to the fea- 
ture in the evolutionary LMC tracks in Fig.|6] see for example the 
50 Mq track for 550km/s. At Galactic metallicity the blueward 
evolution does not appear in the models. Nevertheless, the area 
spanned by the isochrones extends over a wide range of effective 
temperatures. 

Between about 5 and 10 Myr the most massive nonrotating 
stars have evolved off the main sequence. However, at low metal- 
licity, the most massive fast rotators, which undergo chemically 
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Fig. 7. Isochrones for SMC (left panels), LMC (center) and galaxy (right) composition. Each panel shows isochrones for a given age 
and initial rotational velocity ranging from to 540 km/s in steps of 20 km/s, with an increased resolution of lOkm/s between 350 
and 450 km/ s (grey lines). In red we show the isochrone for non rotating models, for other colors see the legend and the main text. 
We show the area spanned by the zero-age isochrones as a light blue band. The age in the second row has been chosen such that the 
isochrones span the largest range in effective temperature, 1 Myr earlier is shown in the first row, 1 and 5 Myr later are shown in the 
third and fourth row. 



homogeneous evolution are still in their main-sequence phase at 
this time, forming a blue straggler-like blue population (see the 
bottom panels of Figs|7]i. If homogeneously evolving stars exist 
in nature they would most likely be found in low metallicity star 
clusters with ages between 5 and 10 Myr. 

4. Surface Abundances 

4.1. Abundances as a function of time 

A direct observable consequence of rotationally induced mixing 
is the enrichment or depletion of certain elements in the atmo- 
spheres of main-sequence stars. In Fig.l9]fmwe show the evolu- 
tion of the surface abundance of various elements as a function 
time. The effects of rotational mixing are more pronounced at 



lower metallicity, at higher masses and for higher rotational ve- 
locity (see also Fig.|4]and Sec.[3]l. 

As is usual in observational work, we express the surface 
abundances relative to the abundance of hydrogen. When stars 
become significantly hydrogen depleted at the surface, using 
hydrogen as a reference element may not be the most logical 
choice. Changes in the abundance may partially reflect changes 
in the reference element hydrogen. We plot the abundances in 
red when these effects become important (i.e. when the helium 
mass fraction at the surface becomes larger than 40%) . 

Most of our stellar models evolve to the red supergiant stage 
directly after the end of core hydrogen burning. This leads to a 
large vertical step in the surface abundances of many elements 
in Fig.l9lfm which is due to the convective dredge-up in the red 
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supergiant stage. In the following, we discuss the evolution of 
the surface abundances of various groups of elements, focusing 
on the changes occurring over the course of the main-sequence 
evolution. 



4.1.1. Helium 

Even though helium is the main product of hydrogen burn- 
ing, the abundance of helium at the surface remains remarkably 
constant in most evolutionary tracks during the main-sequence 
phase. For the 12 models the enhancement is less than 0.2 
dex. Only the fast rotators of 30 and 60 Mq show significant he- 
lium surface enrichments, especially at low metallicity (Fig. |9]- 
fTTI ). These behaviors can be understood as the combination of 
two effects. Firstly, the production of helium occurs on the nu- 
clear timescale. This occurs slower than, for example, the pro- 
duction of nitrogen or destruction of Li. Secondly, with the pro- 
duction of helium a steep gradient in mean molecular weight 
is established at the boundary between the core and the enve- 
lope. Such gradients inhibit the efficiency of mixing processes 
and prevent the transport of helium to the surface. Significant 
amounts of helium can be transported to the surface in models 
were mixing processes are efficient enough to prevent the build- 
up of a chemical barrier between core and envelope. 



4.1.2. The fragile elements, Li, Be, B and F 

The elements Li, Be, B are destroyed by proton captures at tem- 
peratures higher than 2.5 x lO^' K for lithium, 3.5 X 10^ K for 
berylHum, 5 x lO^K for boron (iMcWilham & RauchlllOOl p. 
123). These elements can only survive in the outermost layers. In 
rotating stars the surface abundances of these elements decrease 
gradually as the outer layers are mixed with deeper layers in 
which these elements have been destroyed. The decrease is most 
rapid for the most fragile element, lithium. The surface abun- 
dance of fluorine behaves similarly. It can survive temperatures 
of up to about 20 x 10''K. We note that the rates o f the reactions 
in wh ich fluorine is involved are quite uncertain (lArnould et al.l 
[T999I) . 

In addition, stellar winds can affect the surface abundances 
especially for the most massive stars at high metallicity. Due to 
mass loss deeper layers are revealed in which the fragile ele- 
ments have been destroyed. For nonrotating stars a sudden drop 
in the surface abundances of the fragile elements is visible in 
Fig.l9lfm when mass loss has removed the layers in which these 
elements can survive. In rotating models the change is more 
gradual. 

After the end of the main-sequence, dredge-up leads to a 
further decrease of the surface abundances of these elements. 
However, for some models with masses between 20 and 40 
we find that small amounts of Li and Be are produced in the hy- 
drogen burning shell by the p-p chain. As a result of the interplay 
between a convective zone on top of the shell source and the con- 
vective envelope reaching into these layers after hydrogen ex- 
haustion, we find that the surface abundances of these elements 
can be significantly increased in these models. Even though the 
possible production of these elements by massive stars has sev- 
eral interesting applications, the robustness of the predictions for 
these elements requires further investigation. 



4.1.3. Carbon, nitrogen, oxygen and sodium 

Carbon, nitrogen and oxygen take part as catalysts in the con- 
version of hydrogen into helium. Although their sum remains 
roughly constant, nitrogen is produced when the cycle settles 
into equilibrium, at the expense of carbon and, somewhat later, 
oxygen. The CN-equilibrium is achieved very early in the evo- 
lution, before a strong mean molecular weight gradient has been 
established between the core and the envelope. Therefore, ro- 
tational mixing can transport nitrogen throughout the envelope 
to the stellar surface. While the surface abundance of nitro- 
gen increases, carbon is depleted. Since full CNO-equilibrium 
is achieved only after a significant amount of hydrogen has been 
burnt in the core, changes in the oxygen surface abundance are 
only found in later stages in the more massive and faster rotating 
models. 

Sodium is produced in the extension of the CNO-cycle, the 
NeNa-chain. The changes in the surface abundances of sodium 
resemble the changes is nitrogen, even though the enhancements 
are smaller and appear a little bit later at the stellar surface, see 
Fig.mi] 

The relative increase of the nitrogen abundance depends on 
the initial amount of carbon available for conversion into nitro- 
gen. These amounts are different for the different mixtures. The 
C/N ratio in the SMC and LMC composition are similar (7.4 and 
7.1, respectively), while the ratio in our Galactic composition is 
smaller (3.1). In Fig.|8]we show the evolution of the surface ni- 
trogen abundance in models of 15 and 40 M© with initial rotation 
rates of and 270km/s for SMC, LMC and Galactic composi- 
tion. The initial abundance of nitrogen in the Magellanic clouds 
is much lower than in the Galactic mixture. Nevertheless, the ro- 
tating SMC and LMC models reach surface nitrogen abundances 
at the TAMS which can be higher than in the initial Galactic mix- 
ture. The steep rise at the end of the curves is due to the dredge- 
up when the star becomes a red supergiant. The TAMS is located 
at the base of the steep rise. 
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Fig. 8. Nitrogen as a function of time at SMC, LMC and Galactic 
metallicity. The models are of 15 and 40 M©. Full lines repre- 
sent nonrotating models, dashed Unes models rotating initially 
at ~270km/s. 
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4.2. Surface abundances for selected ages 

Due to the effects of rotation, the age or mass of a star can 
no longer be determined uniquely from its position in the 
Hertzsprung-Russell diagram by fitting evolutionary tracks or 
isochrones. This holds especially for stars that are located near 
the zero-age main sequence. A young, unevolved, slowly rotat- 
ing star may have the same effective temperature and luminosity 
as a less massive, fast rotating, evolved star. A determination of 
the surface composition and to some extent the projected rota- 
tional velocity may help to lift this degeneracy. 

In the panels of Fig. n"2lfT4l we show which abundance dis- 
tributions are predicted at a given time, as function of the con- 
sidered mass and rotational velocity. Models of more massive 
stars show more pronounced surfaces abundances changes for 
a given age. The kinks in the lines of 400km/s at LMC and 
Galactic composition occur due to difficulties in the interpola- 
tion between stars that follow blueward and redward evolution- 
ary tracks effects. 

These plots also show the effect of rotation on the main se- 
quence life time. For example. Fig. [14] shows that the lines rep- 
resenting models with initial rotation rates of 400 and 500 km/ s 
predict that stars of 50 and 60 are still present at 4Myr, 
whereas nonrotating models predict that they have left the main 
sequence. On a much smaller scale this effect is also visible if 
one compares lines for rotation rates between and 300 km /s 
based on models that follow normal evolutionary tracks. 

The steep drop in the boron abundance along the isochrones 
based on nonrotating models is related to the stellar wind mass 
loss. For example, at LMC composition, the drop occurs around 
35 Mq, indicating that the winds efficiently remove the outer lay- 
ers of stars of 35 Mq and higher (see also Sec. l4.1.2] i. 

5. Summary 

In this paper we have presented an extensive grid of models for 
rotating massive main sequence stars. We provide three sets of 
initial compositions that are suitable for comparison with early 
OB stars in the SMC, the LMC and the Galaxy. The models 
cover the main-sequence evolution and in most cases these have 
been computed up to helium ignition. In terms of overshoot- 
ing and rotational mixing efficiency, our models have been cal- 
ibrated against the FLAMES survey of massive stars. We are 
using a new method to calibrate the amount of overshooting 
that makes use of the observed drop in projected rotation rates 
for stars with surface gravities lower than \ogg - 3.2 dex. 
Interpreting this drop as the end of the main sequence, we find 
an overshooting parameter of 0.34 + 0.1 pressure scale heights. 

We have also presented a detailed set of isochrones based 
on models of rotating massive main sequence stars. Whereas 
classical isochrones can be represented by a single line in the 
Hertzsprung-Russell diagram, the isochrones of rotating stars 
span a wide range of effective temperatures at a given luminosity. 
Therefore, the mass and age of an observed star can no longer 
be determined uniquely from its location in the Hertzsprung- 
Russell diagram. We also provided detailed predictions for the 
changing surface abundances of rotating massive main sequence 
stars. While we believe that the data provided here can be use- 
ful to ma ny future studies of massive stars, we make use of it in 
(iPaper Ilh for undertaking a quantitative test of rotational mixing 
of massive stars in the LMC. 
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Fig. 9. Change of the surface abundances of helium, lithium, beryllium and boron along evolutionary tracks of SMC composition 
as a function of time, expressed as a fraction of the main-sequence lifetime. Different line styles represent different initial rotational 
velocities (see legend). For models where the surface helium mass fraction exceeded 40%, the abundances are plotted in red, to 
indicate that the reference element hydrogen is depleted significantly (see also Sec.|4ll. 
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Fig. 9. Surface abundances of rotating models at SMC metallicity, continued. Shown are carbon, oxygen, nitrogen, fluorine and 
sodium (from top to bottom) as a function of time, expressed as a fraction of the main-sequence lifetime. 
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Fig. 10. Like Fig.|9l for stellar models of LMC composition. 
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Fig. 10. continued; Like Fig.|9] for stellar models of LMC composition. 
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Fig. 11. Like Fig.|9]for stellar models with a Galactic initial composition. 
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Fig. 11. continued: Like Fig.|9] for stellar models of Galactic composition. 
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Fig. 13. Surface abundances for a fixed age, as a function of initial mass and rotational velocity, for LMC metallicity. From top to 
bottom are shown helium, boron, carbon, nitrogen, oxygen and sodium abundances. 
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I. Brott et al.: Grids of Evolutionary Models and Isochrones 
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Fig. 14. Surface abundances for a fixed age, as a function of initial mass and rotational velocity, for Galactic metallicity. From top 
to bottom are shown helium, boron, carbon, nitrogen, oxygen and sodium abundances. 



21 



